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ABSTRACT 

We present results of an analysis of the J, H, and K s 2MASS images of 139 spiral 
edge-on galaxies selected from the Revised Flat Galaxies Catalog. The basic structural 
parameters scalelength (h), scaleheight (zq), and central surface brightness of the stellar 
disks are determined for all selected galaxies in the NIR bands. The mean relative 
ratios of the scaleheights of the thin stellar disks in the J:H:K S bands are 1.16:1.08:1.00, 
respectively. Comparing the scaleheights obtained from the NIR bands for the same 
objects, we estimate the scaleheights of the thin stellar disks corrected for the internal 
extinction. We find that the extinction-corrected scaleheight is, on average, 11% smaller 
than that in the K-band. Using the extinction-corrected structural parameters, we find 
that the dark-to-luminous mass ratio is, on average, 1.3 for the galaxies in our sample 
within the framework of a simplified galactic model. The relative thicknesses of the 
stellar disks zq/H correlates with their face-on central surface brightnesses obtained 
from the 2MASS images. We also find that the scaleheight of the stellar disks shows no 
systematic growth with radius in most of our galaxies. 

Subject headings: galaxies: structure, galaxies: spiral, galaxies: photometry, galaxies: 
halos, infrared: galaxies 



1. Introduction 

Study of the parameters of galactic components via decomposition of the rotation curves is 
often ambiguous in the case of spiral galaxies. The same rotation curve often may be explained by 
models with very different parameters of the components. The dark-to-light mass ratio may vary 
several times between the models explai ning well the same rotation curve, especially if the latter is 



smoo th and featureless (see for example Ide Blok et al. 



2001 



Bizyaev fc Zasov 



2002 



Dutton et al 



2005). Incorporating information about the stellar velocit y dispersion from observations helps to 



significantly constrain the masses of galactic components (jBottema 



1993 



Khoperskov et al. 



2001 



and relieve some ambiguities that appear in the modeling of rotation curves, but observing the 
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distributions of the velocity dispersion is a rather diffic ult task. The stell ar disk thickness can be 
used as an alternative to the stellar velocity dispersion (jZasov et al.ll2002l ). Studies of the vertical 
structure of stellar disks alone allow the constraint of some dark h alo parameters (|Bahcaillll984l ; 
Zasov et al.lll99ll . 120021 ; iBizvaev et al.ll2003l ; iBizvaev Kaisinll2004l ) . Note that not all galaxies are 
suitable for such studies: inferring the structural parameters for a spiral galaxy with a significant 
bulge from photometric data alone is dangerous and should be supplemente d with concerns of 
different mass-to-light ratio for different galactic subsystems (see discussion by lAbadi et al.ll2003l . 
regarding difference between photometrically and dynamically identified subsystems). 

Infrared observations are crucial for studies of the structure of edge-on galaxies. The dust 
attenuation makes the extraction of the structural parameters difficult from the optical data. In 
the NIR bands the galaxies are much more transparent, whi ch facilitates studies of the structure of 
galaxies. The all-sky 2MASS survey (ISkrutskie et al.ll2006l ) offers a good opportunity to increase 
the number of edge-on galaxies available for studies in the near-infrared bands. Although faint 
parts of the galaxies are not seen in the 2MASS images, the thin disks of the galaxies are obtained 
with sufficient signal-to-noise ratio (S/N hereafter). 

In this paper we focus on the determination of stellar disk structural parameters corrected for 
the internal extinction. One approach to the problem is to perform the modeling of the r adiation 
trans fer in 3-D stellar-gaseous disks (like in 



Xilouris E. et al 



1999 



Kylafis et al. 



2001 



Bianchi 



20071 ). Such studies require rather deep images, although ther e are attempts to apply this kind 
of analysis to larger sample of moderate quality (jBizyaevi 120071 ). In the present paper we develop 
another approach to study the extinction-corrected structure of the stellar galactic disks comparing 
the parameters of the disks obtained from different NIR bands. 

In Section 2 we discuss the sample selection and evaluation of the stellar disk's structural 
parameters, consider reliability of the parameters, and study the radial variation of the disk scale- 
height. In Sections 3 and 4 we show how the comparison of the parameters obtained in the three 
NIR bands can be used to estimate the extinction-corrected scaleheight in the galactic disks. Sec- 
tion 5 considers the connection between the galactic disk thickness with parameters of the dark 
halos around the spiral galaxies in the frames of a simplified model. Section 6 summarizes results 
of the paper. 



2. Sample of Galaxies and Analysis of the Images 



As a pre-defined sample of edge-on galaxies, the Revised Flat Galaxies Catalog (RFGC here- 
after, iKarachentsev et al.lll999l ) is chosen. The catalog contains visually-classified extragalactic 
objects with major-to-minor axes ratio a/b > 7. To secure enough pixels for studying the vertical 
structure, we select objects with the major axis greater than 1 arcmin. This size is estimated 
directly from the 2MASS images at the level of S/N ~ 3. Note that this size is significantly less 
than commonly used diameter D25. As a result, we select 139 spiral edge-on galaxies that have 
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images in all three 2MASS bands: H, J, and K s 



We apply the technique described by iBizvaev Mitronoval (|2002l ) to obtain the structural 



parameters of stellar disks of the galaxies disks. This method is based on the analysis of photometric 
profiles drawn parallel to the major and minor axes of galaxy at a one-pixel interval. The volume 
brightness in the stellar disks is assumed to change in the radial r or vertical z directions as follows: 

PL(r, z) = p L0 exp(-r/h) sech 2 (z/z ), (1) 

where h and zq are the scalelength and scaleheight of the disk, respectively, and Iq = J p(0, z)dz is 
the central surface luminosity of the face-on disk, which corresponds to the central surface brightness 
^o- The model photometric profiles are obtained by the integration of equation (p} along the line- 
of-sight and then by convolution with the instrumental profile. The best-fit parameters h, zq, and 
Iq are estimated for each radial or vertical profile, and then their median values are taken as the 
resulting scales and central surface brightness of the stellar disks, whereas the standard deviations 
of h, zq, and Jo estimated from different profiles are taken as uncertainty of the parameters. The 
central regions of the galaxies in which bulges may exist are excluded from our analysis. The 
parameters h, zq, and pq are obtained in the J, H, and K s bands for selected 139 galaxies. The 
K s -based parameters are shown in Table 1. 

Figure 1 shows histograms of distribution of the scale heights (top panel) and scale lengths 
(bottom panel) estimated for our program galaxies. In both panels, the J- and H-band scales are 
normalized by those in the K s band. It can be noticed that the scaleheights look systematically 
larger, on average, in the J band than those in the K s for the same galaxy. We will discuss this 
effect below: in §[3] using our observational data, and in § 12.21 with the help of artificial data. 



2.1. Radial Variation of the Disk Scale Height 



As in IBizvaev Sz Mitronoval (120021 ) , we investigate how the scaleheight zq of the stellar disks 
in our galaxies change with the distance to the center r. Our approach to the estimation of the 
zq is independent of the apriori assumptions about behavior of the scaleheight with radius, in a 
contrast with 3-D modeling of the radiation transfer in the stellar and dust disks. Since we analyze 
the vertical profiles separately from each other, this approach is insensitive to possible bending of 
the galactic disks. Note that if the resulting scaleheight zq(t) demonstrates a systematic growth 
with distance to the center, it suggests that the real radial gradient of the disk scaleheight may be 
a few times greater due to the projection effects. 

We fit Zq(t) with a linear function of radius r and find the radial gradients of zq (expressed in 
units of Zq) per scalelength, see Figure 2. Different lines in Figure 2 show histograms of distribution 
of the gradient in the J, H, and K s bands in our whole sample. The median gradient {(Izq/ dr)-(h/ zq) 
in the K s is essentially zero (-0.01), as well as in the H (0.) and J (0.01) bands. 



Note that some galaxies in histogram in Figure 2 indeed show noticeable radial gradients of the 
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scaleheight. In combination with projection effects, it suggests significant gradients {dzQ / dr) ■ {h / zq) 
in a small fraction of our galaxies. 



2.2. Reliability of the Structural Parameters 

The dust disks embedded into the stellar components of the galaxies may significantly distort 
the structural parameters of the stellar disks estimated even from the NIR images. We check 
reliability of the estimated parameters via the following simulations. A set of artificial images is 
created using equation ([T]) and typical parameters h and zq. The surface brightness /iq is replaced 
by the luminosity surface density Iq in our simulations. Each stellar disk contain s an embedde d 



dust disk inside with scalelength of 1.5h and scaleheight of 1/3 zq (according to iBianchil 120071 ). 



The face-on opacity of the dust disk corresponds to Ay = 1 mag. This value should be close to 



the upper limit for t he face-on extinction in spiral galaxies (jKuchinski et alJll99a : IXilouris E. et al 



1999 



Bianchil 120071 ) , Finally, the model disks are randomly inclined by 85-90 degrees and their 
images (i.e. 2-D projections to the plane of sky) are obtained via integration along the line-of-sight. 
The images are then convolved with the instrumental profile typical for the observations. Then 
the poisson and gaussian noise patterns are added to create sets of typically thirty artificial images 
with the same parameters but with different random noise pattern. The images are analyzed by the 
same way as the real galaxies; the mean values of h, zq, and Iq and their r.m.s. are estimated and 
then compared with the input parameters. Figures 3 and 4 give examples how well the structural 
parameters can be extracted from the artificial edge-on disks made with stars and dust should they 
be observed in the J and K s NIR bands. Although the scalelength and the central surface brightness 
can be recovered with a fairly good accuracy from both J and K s bands, the scaleheight in the J 
band looks systematically thicker due to the dust absorption effects only. The K s scaleheight in 
this case is found to be about 15% greater than the model input scaleheight. 

We check how the chosen values of the model parameters affect our conclusions. The amount 
of the dust opacity is responsible for increasing the estimated scale height and makes it thicker by 
10-15% in the K s if Ay is of the order of 1 mag. Much higher values of the central dust opacity, 
more than 3 mag, make stellar disks look about twice as thick even in the K s band. At the same 
time, the disk scalelength shows relatively small variations in this case. Variations of the dust 
scalelength between h and 2h do not change the estimated structural parameters significantly. The 
dust disk thickness affect the estimated stellar disk thickness in the most power when it is about 
a half of zq. Once the dust disk is getting comparable in thickness with the stellar one, it affects 
mostly the stellar disk surface brightness, and not zq. 

The inclination of the stellar disk to the line-of-sight is an important parameter and it increases 
the stellar disk thickness estimated in our modeling by more than 15%, if the inclination is getting 
to be greater than 5 degrees. On the other hand, the relatively small inclination has almost no 
effect on the recovered stellar disk scalelength and central surface brightness. More than 5 degrees 



inclination in galactic disks would create problem with estimation of zq (see also Ide Grijs et al 
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19971 ). but such galaxies should fail the RFGC catalog criteria and we don't expect to have many 
of highly inclined disks in our sample. If the size of the modeling disk is so small that its true 
scaleheight is much less than the seeing, the resulting zq gets independent of the original disk 
scaleheight. Our selection of objects for the analysis by size of their major axis prevents such 
situation. Given all other parameters fixed, the recovering scalelength, scaleheight, and surface 
brightness are independent of the numerical value of h, zq, and Iq unless zq is getting smaller than 
a pixel in the model images. 

In general, our modeling suggests that we can estimate the structural parameters with of order 
of ten percent precision using our method. At the same time, the sensitivity to the variations 
of the scaleheight within the same galaxy with different amount of extinction (like in the case of 
observations in different NIR bands) should be by the order of magnitude better. 



The Basic Structural Parameters in Different NIR Bands 



As we can see from our observational results shown in Figure 1, the stellar disks look thicker 
in the J band and thinner in K s for the same object, in agreement with our modeling from § 12.21 
Two factors may be responsible for such a difference in real galaxies: the vertical gradient of 
stellar population properties and the dust attenuation. Let's consider the effects of the radial and 
vertical gradients of stellar population in galactic disks. The radial grad ients of metal l icity in the 



thin disks are of the ord er of -0.07 dex/kpc for the late-type galaxies (Ivan Zee et al.lll99S), an d 



even less in our Galaxy (jRolleston et al.l I200Q ; iDaflon Cunha 



2004 



Allende Prieto et al 



20061 ). 



According to the population synthesis models by lBruzual fc Charlotl ([2003), this creates gradients 
in (J-K s ) color of the order of -0.015 mag/kpc (i.e. the periphery is bluer). Hence, the scalelength 
in J should appear 4 % longer than that in K for a typical scalelength values of 3 kpc. A radial 
age gradient in combination with the metallicity change may rise this scalelength change to 5%. 
As we see, these numbers are within our accuracy of estimation of the scalelength. The vertical 
gradients of metallicity and age, and hence contribution of the stellar population ef f ects to the 



verti c al color variations within t he thin stellar disk, are usually small (|Tadrossl 120031 ; ISeth et al 
20051 ; lAllende Prieto et al.l 120061) . Even if we assume the large values for the vertical gradient 



of [Fe/H] as in Ide Griis Peletierl (|2000l ). -0.2 dex/kpc, this would create less than 1% longer 
scaleheights in the J than that in the K s due to the stellar population gradients. The amplitude of 
the changes of scalelengths and scaleheights between the H and K s bands is roughly a half of that 
between the J and K s . 

We conclude that the scaleheights of the pure thin stellar disk (i.e. without dust) esti- 
mated from the J, H, and K s images should appear the same. All significant scaleheight vari- 
ations have to be addr e ssed to the reddening effects. This is in agreement with conclusions by 
Dalcanton &: Bernstein! (120021') . Note that in some cases of extraordinary objects like a superthin 



galaxy UGC 7321 (jMatthews Woodll200ll ). the reddening alone cannot explain the observed op- 



tical color (B-R), (R-H) variations in the vertical direction and vertical gradients of metallicity 
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and age are required for the stellar population. At the same time, iMatthews Woodl (|200ll ) note 
that the superthin's H-K color gradients might be explained by the wavelength-dependent dust 
absorption alone. 



4. Structural Parameters of the Thin Stellar Disks Corrected for Extinction 



Following iKylafis Sz Bahcalll ()1987l ). we consider a vertical luminosity profile of an edge-on 
image of a stellar disk with the volume brightness pl{z) = Plo exp(— z/zq). Embedded co-planar 
dust disk with its scaleheight Zd < z absorbs and scatters the stellar light in dependence of the dust 
ex tinction coeffic ient k\ = K,\Qexp{—z/zci). For simplicity, we neglect the dust scatter appealing 
to iBianchil ([20071 ) conclusion that the structural parameters of stellar disks in 3- D modeling of the 



light distribution with or without scattering are not significantly different. As in lKvlafis Bahcall 



(jl987l ). the brightness distribution along the minor photometric axis (z) in such a disk is I{z) 



Iq{z) 1 ~ cx p(~ t ( z )) ; w ]2 ere Jq(z) is the brightness profile in the case of zero extinction, and r is the 
optical depth across the whole edge-on disk, which also depends on the vertical distance over the 
galactic plane z. For small r, I{z) ~ Jq(1 — 0.5 t(z)). The observed scaleheight z e of the stellar 
disk with added extinction is connected with its projected brightness along the selected profile as 
I(z) = I ce exp(— z I z e ) . Given Iq{z) = I c exp(—z/z s ), we find that l/z e ks 1/z s (1 — t(0)). Therefore, 
estimating the thickness of the stellar disk with dust from its edge-on profiles we expect to find 

z eX = z s (l + C x r v ), (2) 

where z e \ is the z e estimated from the disk's brightness profile at a certain photometric band, and 
the optical depth in the plane of the disk in certain photometric band is t\ = C\T y. We assume 



that Ca equals to 0.276, 0.176, and 0.112 in the J,H, and K s bands, respectively (ISchlegel et al 



19981 ) . Since the V-band optical depth Ty and z s are fixed for the same disk, we use linear regression 



and fit a line to three data points (for J,H, and K s values) at C\ - z e \ diagram to find z s . 

For the convenience of comparison of the vertical scales, we normalize the J and H scaleheights 
in each galaxy by their K s -band scaleheight. The upper panel in Figure 5 compares the normalized 
scaleheights. Each line in the panel connects three thicknesses for one galaxy. The tendency to 
observe the stellar disks thicker in the J than in the K s is well seen. The histogram in the lower 
panel in Figure 5 shows the distributions of all normalized zq in three NIR colors. The vertical 
thick line designates all scaleheights in the K s because of the normalization. The mean ratio of Zq 
in the J and H to that in K s is 1.16 and 1.08, respectively. 

We calculate the scaleheight z s corrected for the extinction (or extinction-free hereafter) using 
equation (|2|), which gives z s < Zq[K s ). The distribution of z s normalized by the zq(K s ) is also 
shown in the lower panel in Figure 5. The median z s /zo(K s ) is about 0.89, i.e. the thin stellar 
disks in our galaxies are 11% thinner than it can be estimated from the analysis of the vertical 
profiles from the K s -band images. The extinction- free scaleheights are given in Table 1. 
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We observe no systematic variations in the scalelength between the NIR bands, as it can be 
seen in Figure 6, which is the same as Figure 5 but drawn for the scalelengths. This is not a surprise 
because our simulations shown in Figures 3 and 4 suggest rather small systematic variations of the 
radial scalelength among the NIR bands, and proximity of the h estimated from the J, H, and K s 
bands to its extinction-corrected value. Using the scalelength h from the K s band and extinction 
corrected scaleheights we find that the mean radial-to- vertical scale ratio in the thin extinction-free 
stellar disks of galaxies in our sample is 5.6. 



Thickness of the Stellar Disks and Mass of the Dark Halos 



The thickness of the stellar disk is sensiti ye to the gray itational potential of the disk + bulge 
+ halo system in the vertical direction (e.g. iBahcalll Il984l ) . Therefore, the relative thickness of 
the stellar disk zq/H may provide some information about the spherical-to-disk mass ratio. As a 
first step, we consider a very simplified model of a galaxy that consists of a stellar exponential disk 
embedded into a spherical dark halo. 



The total mass of the exponential stellar disk with the central surface density T,q is 

M d = 2vrS /i 2 . 
The total mass of all components of the galaxy is 

M t ~ V 2 h, 



(3) 



(4) 



where V is the disk's rotational velocity. Following iZasov et al.l (|199ll ) and iKregel et al.l (|2005l ). 
we assume that the disk thickness depends on the local surface density £ and vertical velocity 
dispersion a z : 

o*/GZ (5) 



20 



and the ratio of the vertic al-to-radial velocity dispersions i s a z /a r = const. This is a rather coars e 



but fair appr o ximat ion, see lPolvachnko fc Shukhmanl (jl977l ); lBottemal (|l993l ): lGerssen et al.l (|2000l ) 
Kregel et al.1 ifcoOfJ ). 



We incorporate an assumption of the marginal stability of the stellar disk, i.e. 

o r = Q ■ 3.36S/a3, 



(6) 



where the Toomre parameter Q = const (jToom 



and ae ~ V/h for the flat rotation curves (see 
and £ from equations ([3D, (UD, and 
~ (M d /M t )h. Therefore, 



196 4)- Here se designate s the epicyclic frequency, 



Binney Sz Tremainelll987l . p. 121). Substituting a z 
into ©, we find z ~ Q 2 S 2 /(a3 2 GS) ~ Zh 2 /(V 2 G) 



z /h ~ (M d /M t ). 



(7) 
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Numerical N-body simulations (jZasov et al.lll99ll ; iMikhailova et al.l l2000) in which the spher- 
ical s ubsystem was int r oduce d as a fixed potential confirm the equation ©. Figure 7 is adopted 
from IMikhailova et al.l (|2000l ) with additions from later simulations by Khoperskov (private com- 
munication) . Figure 7 shows that the stellar disk thickness Zo/h is sensitive to the relative mass of 
the spherical component M s normalized by the Both M s and are determined within the 
limits of the stellar disk (four scalelengths in the simulations). Circles in Figure 7 correspond to 
different numerical models with various parameters of the disk and dark halo subsystems, as well 
as with different bulges. The upper curve corresponds to the bulgeless models, whereas the lower 
one is for the models with some bulge contribution. 

Contemporary and more realistic N-body simulations include "live" halos that respond to the 
gravitational potential of disks. The relation zo/h versus M s /Md should look very similar in those 
simulations as in Figure 7 for the galaxies in which the spherical subsystems dominate by mass 
(i.e. for large M s /Md). For the spherical subsystems with masses M s < Md, Figure 7 illustrates 
an extreme case, and its top left part (large z$/h and small M s /Md) demonstrates the lower limit 
for mass of the spherical subsystem given disk thickness Zo/h. Preliminary results from N-body 
simulations (Khoperskov, private communication) with live halo suggest not essential corrections 
to M s /Md inferred from Zo/h according to Figure 7 (of the order of 15%). We encourage reader 
to assume our further estimations of halo-to-disk mass ratio made from the disk thickness as lower 
limits for the real M s /M^. 

Following Figure 7, we calculate the distribution of the dark halo mass Mhalo expressed in the 
units of the stellar disk mass for our galaxies. The RFGC members are "flat" galaxies that are 
not expected to harbor large bulges, so we assume here that M^aio = Ms- We find the following 
relation for bulgeless galaxies from Figure 7: M^ / /Md = — 0.9419 + 0.3737(/i/zo)- According 
to this formula, the median value of Mhaio/Md is about 1.3 within our sample. The histogram of 
distribution of M^aio/Md is shown in Figure 8. No any correlations between Mh a i jMd and absolute 
magnitude, maximum of rotational velocity, or radial gradient of the scaleheight is found for our 
galaxies. A tendency to observe higher Mh a i /Md in the objects with larger scalelength h is found, 
although this relation shows large uncertainty. 

Following the same basic assumptions, we can include the central surface density So into 
consideration. As it follows from our definition of the disk and total galaxy mass, Mj/Mt ~ 
T,q}i 2 /V 2 h = T,Qh/V 2 . We perform a linear fit between log V and log h for our galaxies with known 
V and found th at h ~ ^ 1 ' 42 - It se ems more correct t o utilize non edge-on galaxies to study such a 
correlation, and iKormendvl (|1990l ) and IWalkerl ( 19991 ) suggest h ~ V 1A and h ~ V 1 ' 5 , respectively, 
for arbitrary inclined disks. If the latter relation is accepted, it gives M^/Mt ~ T,o/V ' 5 , and in 
combination with equation ([7]) we come to the following relation between the disk thickness and 
its central surface density 

V^-So/Ti 1 / 3 . (8) 

The upper panel in Figure 9 shows the relation between h/z s (i.e. extinction-corrected inverse 
thickness of the disks) and the central surface brightness in the K s band /iq(K). Here /J-o(K) was 
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corrected for the extinction according to formula (|2|) . The correction is rather small, and its median 
value for our sample of galaxies is about 0.1 mag/arcsec 2 . 

The lower panel in Figure 9 demonstrates the plot of h/z s versus the central surface density 
estimated from observations as 2.30(M/L)2.5 1 2 24 ~ Mo . The numerical coefficient in the formula 
corresponds to the K s -band solar absolute magnitude 3.33 mag. We assume the constant (M/L#- s ) 
= 1 here. The trend in Figure 9 is in agreement with formula ([8]). The very big scatter of 
points in Figure 9 can be attributed to the relatively low accuracy in the fj,Q, h, and z s , and even 
more uncertainty is due to the non-constant value of the mass-to-light ratio (M/L). The scatter of 
this relation can be lowered dramatically if only large (h > 4 kpc) galaxies with the scaleheight 
gradient in the K s (dzo / zq) / (dr / h) < 0.05 are considered, see Figure 10. In this case the inverse 
disk thickness (h/z s ) depends on (io(K) as 

h/z = 47.557 - 23.069 x + 3.016 x 2 , (9) 

where x = log(2.3 ■ 2.5 1 2 24_w( ^)) = 0.4 * (25.11 — fj,o(K)). The solid curve in Figure 10 represents 
this equation. Apparently, the large spiral galaxies whose stellar disks are not flared can be better 
described by the toy-model developed above. Establishing the correlation between the fio and h/z s 
may be useful for estimating the disk thickness and M^aio/Md in disky galaxies with arbitrary 
inclination. 

Note that our model is oversimplified, and attempts to estimate the Mhaio/^d ratio only. 
Should we need to know the values of Mh a i or taken apart, or structural parameters of the 
dark halo, the rotation curve modeling has to be incorporated. Such additional parameter as the 
disk thickness helps to decrease the ambiguity during the rotation curves decomposition. A special 
attention should be given to correspondence between the photometrically and dynamically inferred 
parameters of galac tic subsystems. An example of N-body simulations of an early-type disk galaxy 



(jAbadi et al.ll2003l ) demonstrates that photometry may fail to trace the dynamically distinctive 
subsystems in the central parts of bulge-harboring galaxies. In our consideration we avoid central 
parts of galaxies, and using the RFGC objects prevents us from considering bulge-dominated stellar 
systems. 

The weakness of the simplified assumptions in the modeling described above encourages us to 
apply a more realistic simulations of the stellar disk vertical structure. Since it requires extensive 
numerical simulations and better quality infrared data, and will gain from the incorporation of 
rotation curves decomposition, we defer such more complex consideration to a further paper. 



6. Summary and Conclusions 

We study how the dust absorption affects the basic structural parameters (disk central surface 
brightness, vertical and radial scales) of stellar disks in spiral galaxies estimated from the NIR 
J,H, and K s images. Using 2MASS data, we compare the structural parameters estimated from 
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different NIR bands focusing on the scaleheight of thin stellar disks. The stellar disks look thinner 
in the 2MASS K s band in the comparison with the H and J. We employ this fact to figure out the 
extinction-corrected scaleheight zq of the thin stellar disk. 

Using 139 relatively large galaxies selected from the 2MASS catalog, we find that the mean 
vertical scaleheight (zo) ratios are 1.16:1.08:1.00:0.89 in the J:H:K s :extinction-free bands, respec- 
tively. This means that the radial scalelength estimated from the K s images is very close to the 
extinction-corrected one, whereas the scaleheight is overestimated by 11% on average. The mean 
extinction correction for the face-on central surface brightness is only about 0.1 mag/arcsec 2 in the 

The median radial-to-vertical scale ratio for our sample is about 6. Using a relation between 
the stellar disk thickness and the halo-to-disk mass ratio obtained from a simplified model, we 
estimate the dark-to- luminous mass ratio in our galaxies within the limits of their optical disks. 
Its mean value is about 1.3. A relation between the stellar disk thickness and its central surface 
brightness is observed for the galaxies from our sample, although it shows significant scatter and is 
affected by uncertainty in additional parameters. This relation is more prominent for large (h > 4 
kpc) galaxies from our sample and can be utilized for coarse estimating the disk thickness and 
relative mass of the dark halo in large spiral galaxies arbitrarily inclined to the line-of-sight. 

We find also that the scaleheight does not change systematically along the radius in most 
of our galaxies. Only a small fraction of galaxies demonstrates a noticeable radial growth of the 
scaleheight. 

This project was partly supported by grant RFBR-07-02-00792. This publication makes use 
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tion. Computational equipment for DB was provided by AAS under a Small Research Grant. We 
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Table 1: The Structural Parameters of the Galaxies 



Name 


h(K) 


dh 


zO(K) 


dz 


Mo(K) 


d^ 


z (cor.) 


zq (cor) / h 


/izo(K,cor) 


M-dark 




kpc 


kpc 


kpc 


kpc 


mag/D" 


mag/D" 


kpc 




mag/D" 




RFGC0095 


0.72 


0.06 


0.31 


0.08 


17.69 


0.26 


0.27 


0.375 


17.61 


0.06 


RFGC0099 


4.35 


0.06 


0.95 


0.17 


17.51 


0.18 


0.83 


0.191 


17.41 


1.02 


RFGC0102 


1.76 


0.03 


0.48 


0.11 


17.49 


0.20 


0.46 


0.260 


17.46 


0.49 


RFGC0124 


8.15 


0.24 


1.36 


0.19 


17.94 


0.18 


1.12 


0.137 


17.81 


1.79 


RFGC0152 


5.33 


0.06 


0.88 


0.14 


17.73 


0.21 


0.76 


0.143 


17.66 


1.67 



The full table is available on-line 

The table shows RFGC name, scalelength h in kpc and its uncertainty in the K s band, scaleheight zq 
in kpc and its uncertainty in the K s , central surface brightness /j,q in the K s , extinction-corrected 
scaleheight zo(cor.) in kpc, extinction-corrected relative thickness of the stellar disk zo(cor.)/h, 
extinction-corrected central surface brightness (io (cor.) in the K s , and dark-to-luminous mass ratio 

Mdark/Mi um . 
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Fig. 1. — The scaleheights (top panel) and scalelenghts (bottom panel) in the J (dashed line) and 
H (dotted line) bands normalized by those from the K s for our 139 galaxies. The K s band scales 
in both panels equal to unit for all objects. 
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Fig. 2.— Distribution of the radial gradients of the scaleheight per one scalelengths 
(dzo/zo) I (dr/h). The solid, dashed, and dash-dotted lines designate the K s , H, and J gradi- 
ents, respectively. The three short vertical lines above the histogram indicate the median values of 
the gradients in the K s , H, and J bands (from left to right), that are essentially equal to zero. 
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Fig. 3. — The structural parameters of the dusty stellar disks estimated from sets of artificial images 
in the J band are shown by the solid lines with squares for different signal-to- noise ratio (S/N in 
the center of the disks is used). The dotted lines designate the r.m.s scatter of the parameters. 
The dash-dotted horizontal lines denote the input (i.e. true) model parameters. The scalelength, 
scaleheight, and central surface luminosity (Iq, in arbitrary units) are shown in the upper, middle, 
and lower panels, respectively. 
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Fig. 4. — The same as in Figure [3] but for the K s band. 
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Fig. 5. — Top panel: all scaleheights zq normalized by that in the K s . Each a galaxy is shown by 
a dashed line, which connects the zo values in the J, H, and K s (the latter equals to unit for all 
objects). The solid line shows the case of equality of all scaleheights. Bottom panel: histograms of 
distribution of the scaleheights normalized by those in the K s band. Dashed line corresponds to the 
J band, dotted - to the H band. The solid line designates the extinction-corrected scaleheights. The 
numbers under the letters show the mean values of the scaleheights in the corresponding bands. 
The thick vertical line designates formally the scaleheights in the K s that are all equal to unity due 
to the normalization. 
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Fig. 6. — Top panel is the same as in Figure 5 but drawn for the disk scalelengths. Bottom panel 
is similar to the bottom panel in Figure 5 and is plotted for the scalelengths. The scalelengths are 
normalized by those from the K s band. Dashed line corresponds to the J band, dotted - to the H 
band. The thick vertical line designates formally the scalelengths in the K s that are all equal to 
unity due to the normalization. 




Fig. 7. — The relative disk thickness Zq/Ji is shown against the spherical-to-disk mass ratio M s /Md- 
Here M s is the mass of the spherical component (halo and bulge, in a general case), and is 
the mass of the disk component. The M s and are estimate d within the limits of the stellar 
disk (4/i). This figure comes from simplified N-body simulations (jMikhailova et al.ll2000l ) and each 
symbol corresponds to certain numerical model. The models were run with a variety of parameters 
for the disk, dark halo, and bulge. The upper curve corresponds to the bulgeless models, whereas 
the lower one is for the models with some bulge contribution. 
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Fig. 8. — The dark halo-to-disk mass ratio for the galaxies from our sample inferred from the 
relative disk thickness. The mass ratio is close to the dark-to-luminous mass ratio since most of 
our galaxies have no prominent bulges. 
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Fig. 9. — Top panel: the inverse disk thickness h/z s corrected for the internal extinction versus the 
central surface brightness in the K s band. Bottom panel: the same relation where hq is replaced 
with the central surface density calculated as 2.30 • 2.5 12 24_Mo , see text. 
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Fig. 10. — The inverse disk thickness h/z s corrected for the extinction versus the central surface 
density of the stellar disks (see text) for the sample of large (h > 4 kpc) and non-flared disk galaxies. 
The solid line shows the best-fit to the data (see text). 



